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1. ORIGIN OF VARIATION 


Variability is a deeply ingrained 
characteristic of nature and mans 
appreciation of variability is the source of 
knowledge, beauty and tolerance. 

Kubis(1966) 


In the search for general trends and laws in Nature, there is a tendency to exorcise 
variation as an artifact, for which imperfection and inaccuracy of observation are held 
resp>onsible. However, there is often a residual variation and unpredictability when mea- 
surement error and external uncontrollable disturbance (noise) have been effectively 
removed. The following principle can be stated: 

Variation is an inherent characteristic of objects and systems — be they 
cells or organisms, cities or countries, rocks or mountains, electric bulbs or 
airplanes — a consequence of the statistical (macroscopically apparently non- 
deterministic) nature of the processes by which the objects and systems are 
formed (or manufactured) from simpler elements. 

Thus, physiological equality among individuals is an impossibility in sexually 
reproducing organisms. Individuals arc “manufactured” based on an elaborate but 
basically random “shuffling and reshuffling of the genetic cards” (the genetic infor- 


ANGELOS C. ECONOMOS and RALPH C. BALLARD • Deparimenl of Biological Sciences, San Jose 
Slate University, San Jose, California 95192, Dr, Economos’s present address is: Laboratoire de Gdnetique, 
Universite Calholique de Louvain, Louvain-Ia-Neuye, Belgium. JAIME MIQUEL • Biomedical 
Research Division, NASA, Ames Research Center, Moffett Field, California 94035. JOHN E. JOHN- 
SON, Jr. • National Institute on Aging, Section on Experimental Morphology, Baltimore City Hospi- 
tals, Baltimore, Maryland 21224; Department of Neurology, Johns Hopkins University School of Medicine, 
Baltimore, Maryland 21205; and Hitachi Scientific Instruments, Rockville, Maryland 20850, 


187 



188 


ANGELOS C. ECONOMOS et at. 


mation in the gametes) drawn from the available “genetic pool” of a population (Dawk- 
ins, 1976). This is the genotypic variation. In highly inbred populations of experimental 
animals, though the genetic design has become relatively fixed, implementation of the 
blueprint in the actual manufacturing of such a highly complex object as an individual 
does not produce identical copies of the prototype. Although variations introduced at 
the various component steps of the manufacturing process may be small, they culminate 
in observable differences among individuals. This is the phenotypic variation. 

2. ANALYSIS OF VARIATION 

Statistical methodology was originally developed to handle measurement error per 
se. Some of the assumptions, however, pertain only to measurement error and not in 
general to the residual variation. In handling another part of variation, i.e., noise, infor- 
mation theory and communication engineering have been very successful — witness, for 
instance, the transmission of information amidst unimaginable levels of noise (com- 
pared to signal magnitude) in planetary explorations. Finally, intensive study of varia- 
tion in technological objects over the past few decades has culminated in the develop- 
ment of a new branch of engineering, reliability engineering (Bazovsky, 1961; Haviland, 
1964; von Alven, 1964; Shooman, 1968; Brook, 1972; Carter, 1972; Jardine, 1973; 
Bompas-Smith, 1973; Gross and Clark, 1975; Kapur and Lamberson, 1977). The aim 
in this discipline is to develop an effective methodology for designing reliable machines 
and other technological objects or systems in spite of variations in characteristics and 
durability of available components. Although in biological and aging research we are 
not expected to design our object of study, use of concepts from reliability engineering 
can be useful. 

Variation in a set of objects or individuals is described quantitatively by the fre- 
quency distribution, i.e., the mathematical function giving the percentage of the popu- 
lation with a certain value of a measured attribute within the range the attribute varies 
over the given group (see Fig. 1 for examples). With the passage of time, wearing-out 
or aging may shift and/or change the “form” of this function. 

Well known is the normal or bell-shaped distribution. Its widespread occurrence 
stems from an interesting property of distributions, described by the “central limit theo- 
rem” (Papoulis, 1965): 



BODY WEIGHT, mg 



Figure I . Frequency distribution of body weights in two inbred populations of young adult animals from 
this laboratory: 235 male fruit flies (32 days old) raised and kept at 2TC and 125 male C57BL/6J mice 
(8 months old) kept at 24 ‘C. 
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organism above a body temperature of roughly 40^*0 due to fast increase with temper- 
ature of the rate of thermal denaturation of protein; higher body temperature, on the 
other hand, increases the thermodynamic efficiency of the body’s biochemical reactions. 
There are also other thermodynamic (c.g., heat transfer), ecological, and evolutionary 
constraints which are in general different in mammals and birds; this may account for 
the differences in the frequency distributions of Fig. 2. 
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Figure 3. Typical examples of asymmetric, log-normal type distributions, of biological and technological 
characteristics, in which an “optimality mechanism” may be operative (see text), (a) Latency of first mating 
in a population of twenty-eight 4-week-old male fruit flies housed individually in vials, measured from the 
moment three 7-day-old virgin females were placed in each male’s vial. Data from Economos et al. (1979a), 
Tables I and II. (b) Area of mitochondria on electron micrographs of cross sections of fibers of the flight 
muscle of 8-day-old houseflies. From Sohal and Bridges (1977), smoothed histogram from Fig. 5. (c) Age 
of mothers in a group that gave birth to 2191 children. Data from Murphy (1954) cited in Calloway and 
Dollevoet (1977), Table 15. (d) Age of Fellows of the Royal Society of London at election. Data from 
Bartlett (1970), Table 2. (e) Diameter of components turned on a hand-operated lathe to a specification of 
0.150 ± 0.001 in., subsequently polished until a specified surface finish was achieved. From Brook (1972), 
smoothed histogram from Fig. 1 1.1. (f) Repair time of a typical (unspecified) piece of equipment. From 
von Alven (1964), smoothed histogram from Fig. 5.10. 
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Figure 2. Frequency distribution of average body temperature in mammals and birds. From Morowitz 
(1968), smoothed histogram from Fig. 5.4. 

If a characteristic is influenced by a large number of mutually indepen- 
dent factors acting additively, no matter what the form of the frequency dis- 
tribution of these factors is, the composite characteristic has a bell-shaped 
distribution. 

Factors often interact, however, or influence a characteristic in nonadditive ways. 
For instance, body weight ( PF) is proportional to the product of three body dimensions 
(£| • Lj ' L^\ its distribution will not be bell shaped (Fig. 1). However, the logarithm 
of body weight is proportional to the sum of the logarithms of the three body 
dimensions: 


PF oc • Lj * £3 log PF oc log £, H- log £3 -f log £3 

By the central limit theorem then, log W will be roughly normally distributed. Very 
often, therefore, the logarithm of a characteristic will be normally distributed; such a 
characteristic is said to be log-normally distributed, a term coined by Gaddum in 1945 
(Oldham, 1968). 

An important property of this distribution is its asymmetry: a large percentage of 
the individuals have a certain value of the studied characteristic, fewer and fewer of 
them having larger or smaller values (but not both). This suggests another mechanism 
by which a log-normal distribution is obtained, operative where there is competition 
among the individuals of a group for a common resource. Many examples of this phe- 
nomenon have been reviewed by Papageorgiou (1976): size of particles formed from 
grinding or breakage processes, bank deposits, size of firms (compared in terms of num- 
ber of employees or annual sales), size of sentences in a text (in terms of number of 
words), size of farms or towns in a country, size of oil reserves, and so on. 

Another general mechanism that leads to an asymmetric distribution similar to a 
log-normal distribution is the following. In some manufacturing processes, a given spec- 
ification of a certain characteristic (often an “optimal” specification by some criterion) 
may be more difficult to be deviated from on one side of the optimal specification within 
the range of permitted values than on the other side of the optimal specification. In the 
examples of Fig. 2, there is the severe constraint of a high bioenergetic cost to the 
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organism above a body temperature of roughly 40‘C due to fast increase with temper- 
ature of the rate of thermal denaluration of protein; higher body temperature, on the 
other hand, increases the thermodynamic efficiency of the body’s biochemical reactions. 
There are also other thermodynamic (c.g., heat transfer), ecological, and evolutionary 
constraints which are in general different in mammals and birds; this may account for 
the differences in the frequency distributions of Fig. 2. 
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Figure 3. Typical examples of asymmetric, log-normal type distributions, of biological and technological 
characteristics, in which an “optimality mechanism” may be operative (sec text), (a) Latency of first mating 
in a population of twenty-eight 4-week-old male fruit flies housed individually in vials, measured from the 
moment three 7-day-old virgin females were placed in each male’s vial. Data from Economos et aJ. (1979a), 
Tables I and II. (b) Area of mitochondria on electron micrographs of cross sections of fibers of the flight 
muscle of 8-day-old houseflies. From Sohal and Bridges (1977), smoothed histogram from Fig, 5. (c) Age 
of mothers in a group that gave birth to 2191 children. Data from Murphy (1954) cited in Calloway and 
Dollevoet (1977), Table 15. (d) Age of Fellows of the Royal Society of London at election. Data from 
Bartlett (1970), Table 2. (e) Diameter of components turned on a hand-operated lathe to a specification of 
0.150 ± 0.001 in., subsequently polished until a specified surface finish was achieved. From Brook (1972), 
smoothed histogram from Fig. 11.1. (Q Repair time of a typical (unspecified) piece of equipment. From 
von Alven (1964), smoothed histogram from Fig. 5.10. 
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Figure 4. Examples of asymmetric log-normal type distributions whose underlying mechanism is not clear, 
(a) Size of lacunae in cross sections of bones. From Yeager ei al. (1975), smoothed histograms from Fig. 
1. (b) Snowfall in New York City in the years 1919-1953. Data from Haviland (1964), Table 10.4. (c) 
Frequency of deaths in a population of 296 male C57BL/6J mice. From Kunstyr and Leuenberger (1975), 
smoothed histogram from Fig. 4. (d) Frequency of deaths in a population of 4627 male houseflies at 25 'C 
ambient temperature. Data from Rockstein and Lierberman (1959), Table ITT. 


Further examples where the above “optimality mechanism” may be operative are 
shown in Fig. 3; this may be less obvious in other examples (Fig, 4), which could also 
be explained by other mechanisms. In some instances the asymmetry of the frequency 
distribution may be exaggerated to the point that no individual has a value of the stud- 
ied characteristic on one side of a certain most common value. Two examples are shown 
in Fig. 5. As is seen by plotting the distributions on a semilogarithmic plot, each can be 
represented by a single exponential function. 


3. VARIATION AND AGING 

An obtrusive characteristic of organisms is that, with the passage of time, they 
suffer progressive deterioration of structure and function leading to senescence and 
death. The frequency distribution of various structural or functional characteristics of 
young individuals is therefore bound to change as they age. One general trend is seen 
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Figure 5. Examples of negalive-exponenlial type frequency 
distributions, (a) Survival time of 256 patients with malig- 
nant melanoma who had metastases at the lime of admission. 
Data from MacDonald (1963) cited in Gross and Clark 
(1975), Table l.I. (b) Rank of citation of authors (number 
of cited papers by an author, single or first authorship only) 
based on 336 references from a monograph. [References 
appearing in a monograph by Papaikonomou (1974) were 
counted and tabulated.] (c) and (d) are semilogarithmic 
plots of (a) and (b), respectively. 
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in a progressive shifting of the most common (or “optimal”) value of a characteristic 
toward less optimal levels. Another trend may be a tendency of variance to increase 
with age. These trends may be apparent in populations in which no appreciable mor- 
tality has yet set in. However, as progressively more weak individuals die, such trends 
may be less clear, because of introduced bias by elimination of the least optimal values 
corresponding to the individuals that have died; thus, under certain circumstances, vari- 
ation may be reduced in old age. Figure 6 illustrates the observed change of the fre- 
quency distribution with age of a structural characteristic of hamsters. Two examples 
of aging changes in the frequency distribution of fine structural characteristics are 
shown in Fig. 7. A striking effect of aging is seen in Figs. 8 and 9. 

A trend seen in the biological example of Fig. 6 is also evident in the technological 
example of Fig. 10. Here, '‘aging” of the capacitors could be considered to be simulated 
by a larger number of dielectric defects (conducting flaws). Such “aging” shifts the 
frequency distribution of the capacitor strength (breakdown voltage) toward lower val- 
ues. Interestingly, however, the form and variance of the distribution remain 
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Figure 6. Frequency distribution of width of glomerular basement membrane of hamster kidneys at three 
different ages: 0, 19, and 70 weeks. Redrawn from McNelly and Dittmer (1976). 


(a) 



(b) 



CIRCUMFERENCE OF 
MITOCHONDRIA, )um 


Figure 7. Effect of age on ultrastructural characteristics, (a) Relative size of mitochondria in cells of per- 
fused livers of young and old C57BL/6J mice. From Wilson and Franks (1975), smoothed histograms from 
Fig. 6. (b) Circumference of mitochondria in human liver cells in two age groups. From Tauchi and Sato 
(1968), smoothed histograms from Fig. 4. 
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Figure 8. Electron micrographs of cross sections of spermatozoa of young (7-day) and old (84-day) flies. Aging has a dramatic effect on the appearance 
and size of the spermatozoa. From Miquel et al. (1972b), Figs. 17 and 18, X 45,000. 
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MITOCHONDRIA "WIDTH," mm 


Figure 9. Frequency distribution of “length” and “width” of the mitochondria of the spermatozoa of young 
and old fruit flies measured manually on the electron micrographs of Fig. 8. (“Length” and “width” were 
measured as shown.) Aging has a striking effect on the distributions, increasing considerably their asym- 
metry and variance. 


unchanged, unlike in the biological example. This discrepancy probably stems from the 
fact that the capacitors were assumed to "age” uniformly (equal increases of the num- 
ber of flaws in all “individuals”), whereas aging may not progress uniformly in aging 
organisms. 

An intriguing contrast with the examples of Figs. 6 and 10 is seen in Fig. 1 1. At 
first sight, the frequency distributions of salaries for assistant, associate, and full pro- 
fessors (presumably young, mature, and middle-age and older individuals, respectively) 
are very similar to the distributions of the glomerular basement membrane width of 
hamster kidneys at corresponding ages (Fig. 6). The similarity extends to both a shifting 
of the distribution and an increase in variance with age. However, unlike in the biolog- 
ical-structure example, where shifting is toward less optimal values, in the example of 
professorial salaries the trend is reversed! Clearly, this state of affairs does not corre- 
spond with aging changes in underlying biological structure and function; it may be an 
example of the effect of memory and learning (experience, information acquisition and 
processing in the brain) in human affairs. 



196 


ANGELOS C. ECONOMOS et af. 




Figure 10. A “simulated-aging” effect on the frequency distribution of breakdown voltage in batches of 
capacitors. Aging is simulated by considering batches of capacitors with increasing numbers of conducting 
flaws in the dielectric. (Top) Cross section of a capacitor. From Haviland (1964), redrawn from Figs. 6.1 
and 6.4. 


4. APPLiCATlONS 

4. 1. Variation in Mouse Liver Cellular and Fine Structure: Effects of Aging. 
Alcohol, and Antioxidants 

Andrew succinctly stated that “ . . . histological changes [in aging] are common 
enough to lead us to anticipate that he who seeks frequently would find changes in 
ultrastructure” (Andrew, 1968). Apparently, this fact has been appreciated, for 
instance, by a number of investigators who have studied liver cell ultrastructure in 
human subjects, mice, and rats over the past two decades (Tauchi and Sato, 1962; 
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Figure 11. Frequency distribution of salaries of professors in Departments of Physiology of American uni- 
versities. From Ganong (1977), smoothed histograms. 


Andrew, 1962; Schmucker et al, 1974; Wilson and Franks, 1975; Fieri et al., 1975; 
Schmuckerand Jones, 1975; Herbener, 1976; Tate and Herbener, 1976). In a previous 
report by our group (Johnson et al., 1978), the effects of aging and antioxidants on the 
fine structure of the mouse hepatocyte were described qualitatively. Here we present a 
quantitative analysis that illustrates application of the concepts discussed above in his- 
tology and aging. (A preliminary summary was reported in Economos et al., 1979b.) 

We studied five groups of male C57BL/6J mice: 8- and 30-month-old controls; a 
group fed an antioxidant-containing diet (“diet 3”), consisting of lab chow containing 
0.07% of tocopherol />-chlorophenoxyacetate and magnesium thiazolidine carboxylate 
from age 26 to 30 months; and two groups given 10% alcohol in their drinking water 
from age 1 2 to 24 months with or without the antioxidants. 

Light microscopic study of liver from the various groups confirmed the finding 
from previous investigations that in normal aging the overall appearance of the liver 
parenchyma does not change dramatically (Andrew, 1962; Tauchi and Sato, 1962), 
apart from a loss of cells reflected in a 20% decrease of liver weight in old controls 
(Johnson et al., 1978). However, alcohol had considerable effects, while the antioxidant 
diet appeared to protect the liver cells to some extent from alcohol. These conclusions 
are supported by Figs. 12 and 13, which show the frequency distributions of cell sizes 
in the five groups (pooled data from four livers from each group, 1(X) cells from each 
liver on four different sections). The distribution of cell sizes in old controls or old ani- 
mals on the antioxidant diet was similar to young controls, with a small tendency toward 
increased cell size that did not affect the average cell size. From observation and mea- 
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Figure 12. Effects of aging and antioxidants (“diet 3”) on the frequency distribution of mouse hepatocyte 
size (largest diameter measured by using a light microscope). 


surement of cell nuclei, it appears that the number of binucleated cells increased slightly 
but not significantly with age from 15.5% in young controls to 16.7% in old controls and 
16.3% in old animals on the antioxidant diet (Table I). This increase was accentuated 
in alcohol-treated animals (22.3% binucleated cells), with some protection by the 
antioxidants (19.0% in the alcohol + diet 3 group). Though alcohol had a very small 
effect on average cell size, the frequency distribution of cell sizes was dramatically 
affected (Fig. 13). Despite the increase in the percentage of binucleated cells, alcohol 
decreased considerably the percentage of larger cells, an effect that was also repressed 
in the group that received antioxidants plus alcohol (Table I). Finally, the average size 
of the nucleus increased slightly with aging, accompanied by a similar increase in the 
average size of the cytoplasm, the nucleocytoplasmic ratio remaining essentially 



Figure 1 3. Effects of alcohol with or without antioxidanU (“diet 3”) on the frequency distribution of mouse 
hepatocyte size (largest diameter measured by using a light microscope). 
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Table 1. Liver Cell Characteristics: Effects of Age, Alcohol, and Antioxidants ("Diet 3") 



Average size (^m) 

% of cells > 26 fim 

% binuclealed cells 

Young 

25.9 

46.2 

15.5 

Old 

26.8 

54.4 

16.7 

Old -1- diet 3 

27.9 

59.6 

16.3 

Old alcohol 

25.9 

46.2 

22.3 

Old -j- alcohol + diet 3 

27.2 

60.4 

19.0 


unchanged, while variance in cell sizes increased, in general agreement with findings in 
the aging human liver (Tauchi and Sato, 1962). 

The electron microscopic study of hepatocytes in the five groups of mice has been 
described elsewhere (Johnson et al, 1978). Figures 14 to 16 show representative cells 
from control and antioxidant-treated mice. Apparently, qualitative changes of hepato- 
cyte fine structure in aging are slight. Alcohol, however, had a dramatic effect on the 
fine structure of liver mitochondria, i.e., mitochondrial clustering and thinning of the 
mitochondrial membrane, were evident. Antioxidants seemed to repress this effect, but 
apparently the cells had fewer mitochondria and considerably more fat droplets. To 
analyze the effects of aging, alcohol, and antioxidants on hepatocyte mitochondria 
quantitatively, we selected at random five electron micrographs from each experimental 
group and measured manually the “length” (I) and “width'" {W) oi each mitochon- 
drion (400-500 total). (When, for irregular mitochondria, a single width could not be 
measured, two widths were measured at about \ length from the mitochondrial central 
point and their average was used.) We then determined a rough but reasonable estimate 
of mitochondrial “area” {A) on the micrograph as the product LY. W, and an estimate 
of one-half the mitochondrial “circumference” (Q the sum L W, Frequency dis- 
tributions of L, W, A, and C were plotted for all groups and are analyzed below based 
on the concepts presented in this chapter. 

Figure 17 shows the frequency distributions of mitochondrial length and width in 
three experimental groups. All distributions had an asymmetric log-normal type form. 
The asymmetry was the smallest for W in young controls; it was somewhat larger for 
W of old controls, a trend abolished by the antioxidant diet. The same relationship is 
apparent for L, but all L distributions had a considerably larger asymmetry than the 
W distributions. This accentuation of asymmetry in the mitochondrial length distribu- 
tion may be an artifact of the measurement procedure. Indeed, the length of a mito- 
chondrion on an electron micrograph is proportional to the product of the organelle’s 
actual length and the angle at which the mitochondrion was cut, unlike for the width. 
This is due to the ellipsoid form of the mitochondrion. (This phenomenon can be better 
visualized by imagining a cut through a cylinder: the width of the cross section of the 
cylinder is equal to its actual width, but the length of the cross section depends also on 
the angle of the cut.) Thus, L, being the product of two distributed characteristics (the 
angle of cutting also has a frequency distribution, because it depends on the random 
orientation of mitochondria in the cell space), has a distribution tending toward an 
asymmetric log-normal type of distribution, although the actual length may have a less 
asymmetric distribution, comparable to that for Ws. 




Figure 14. Electron micrograph of a hepaiocyle from the liver of an 8-month-old control mouse. X 13,650. 




Figure 15. Electron micrograph of a hepatocyte from a 30-monih-old control mouse. Note an apparent 
slight increase in mitochondrial size compared with Fig. 14. X 13,650. 
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Figure 17. Effect of aging and protective effect of antioxidants (“diet 3”) on the frequency distribution of 
mouse hepatocyle mitochondrial “length" (a) and “width” (b). (“Length” and “width” were measured on 
electron micrographs of hepalocytes at 13,650-fold magnification.) 


Following the same reasoning, the estimated mitochondrial area, which is the 
uct of L and W. should have even more asymmetric distributions than either L or IF in 
all groups Figure 1 8a shows this. However, the estimated circumference, being the sum 
rather than the product of L and W, should have a less asymmetric distribution, in 
agreement with the central limit theorem, mentioned above, which dictates that the 
sum of distributed characteristics tends to have a bell-shaped distribution (Fig. 18bh 
Both area and circumference were affected by the antioxidant diet similarly to L an 
W i.e. the youthful distribution patterns were restored. These trends are clearly dis- 
cei-nedin Fig. 19, where the frequency distributions of L.W, LX W. and L + JFare 
plotted for young and old control animals separately. 

Mitochondrial length and area frequency distributions are plotted for four groups 
(young controls, old controls, old -H alcohol, and old -|- alcohol + antioxidants) in Fig. 
20. Paradoxically, alcohol alone changed the distributions only slightly, but when 
antioxidants were given in the diet to protect against the effects of alcohol, the frequency 
distributions were further shifted or their asymmetry toward larger values was accen- 
tuated compared with the old controls and old + alcohol without antioxidants. How- 
ever, as was mentioned above, alcohol had a dramatic qualitative effect on mitochon- 
drial fine structure, probably impairing their function considerably despite a small 
effect on mitochondrial dimensions: the mitochondria became pathologic rather than 


Figure 16 Effect of the antioxidant diet on hcpatocyte fine structure in a 30-month-old mouse; the tendency 
to increased mitochondrial size caused by aging is repressed by the diet. X 13,650. 
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Figure 1 8. Effect of aging and protective effect of antioxidants (“diet 3”) on (a) estimated mitochondrial 
"area” (“length” X “width”) and (b) estimated one-half mitochondrial “circumference” (“length” -H 
"width”). 




40 80 120 160 200 0 40 80 120 160 200 

MITOCHONDRIA "SIZE”, mm or mm^ MITOCHONDRIA "SIZE", mm or mm^ 


Figure 19. Frequency distributions of mitochondrial “length,” “width,” “area,” and “circumference” are 
plotted on the same graph for young and old mice, to compare the various distributions at the same age. 


aged faster. The antioxidant diet may have partially protected the mitochondria against 
this effect of the alcohol, but the cells now had fewer mitochondria (fewer but “health- 
ier”?); the increased fat-droplet volume would make up for the lost mitochondria. 
Therefore, because alcohol appears to damage mitochondrial fine structure (and prob- 
ably function) considerably, the morphometric analysis cannot be applied in this case. 
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Figure 20. Frequency distributions of (a) “length” and (b) “area” of mouse hepalocyle mitochondria in 
four of the five experimental groups, plotted on the same graphs to show the effect of alcohol. 


The qualitative electron microscopic evidence and the quantitative analysis of light 
microscopic information that showed virtual elimination of larger cells by the alcohol, 
an effect repressed by the antioxidants (Fig. 13), indicate that the antioxidant diet pro- 
tected the mitochondria to some extent from the damaging effect of high levels of alco- 
hol in mice. 


4 , 2 . Variation in Vitality and Mortality 

The individuals of a cohort die at different ages, even in inbred strains of animals. 
Typical frequency distributions of ages at death were shown in Fig. 4c and d. Previous 
theories of mortality have generally assumed that the individuals of a population are 
identical and have sought to explain the variation in age at death based on various 
assumptions on the statistical properties of the causes of death. However, variation in 
constitutional and functional characteristics of individuals is a reality that should not 
be suppressed. Even more, it is shown elsewhere (Economos and Miquel, 1979) that 
variation in mortality and physiological characteristics are closely related. The model 
is summarized below. 

A physiological variable measured over the individuals of a population has a fre- 
quency distribution typified by latency of first mating of flies (Fig. 3a). A number of 
physiological functions can be chosen that are important for survival to a given extent. 
A weighted sum of these variables in accordance with their relative importance will 
have a similar frequency distribution. This sum, called vitality, may be considered to 
reflect an individual organism’s capacity to resist death. Because most physiological 
functions decline uniformly with age (Shock, 1960), so will the vitality of each individ- 
ual, the form of its frequency distribution remaining roughly unchanged with age. 
Therefore, the distribution of ages at death can be derived by a linear transformation 
of the distribution of vitalities (Fig. 21). If the distribution of vitality is changed, e.g., 
by eliminating weak or unfit individuals through a physical examination, the distribu- 
tion of ages at death and mortality kinetics will change accordingly (Fig. 22). 
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Figure 21. Hypothetical frequency distribution of vitality in a pop- 
ulation and its relationship with population mortality kinetics. For 
the mathematical model it is assumed that the distribution stands 
still as the population ages, while a life-death threshold increases at 
a rate b (actually the reverse is true). At each age total mortality 
will be proportional to the area under the vitality curve below the 
life-death threshold. (For a detailed description of this model, sec 
Economos and Miquel, 1979.) 



Figure 22. Schematic representation of 
changes in vitality distribution caused by 
selection through a physical examination 
of the “healthier,” more fit individuals 
(medically insured persons and pilots), and 
predicted effects on mortality kinetics, 
which agree with observations. From Econ- 
omos and Miquel (1979). [Cumulative 
mortality is plotted semilogarithmically, as 
proposed by Economos (1979a).] 


(a) (b) 



4.3. Time-Condensing in Experimental Aging Research through the Study of 
Variation 

An important shortcoming of the classical actuarial methodology used in applied 
aging research, where the effect of environmental factors and dietary chemicals (aging 
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rate modulators) is assessed in terms of effects on mortality kinetics, is well known: the 
exp)erimental animals of choice (mice and rats) live too long. It is therefore imperative 
to have a methodology that will enable us to predict the effects on life span in a much 
shorter time than the 2\ and 4 yr that are the average and maximal life spans of mice 
and rats. This is the meaning of “time-condensed studies” (Economos and Miquel, 
1977). Two approaches are conceivable: 

1. Evaluate the effect of an aging rate modulator on the frequency distribution of 
a vitality variable of an experimental population by comparing with a genetically sim- 
ilar control population over an appropriate period of time (a fraction of the life span). 
From this comparison, an estimate for the effect of the aging rate modulator on aging 
and life span can be derived. 

This approach is based on the hypothesis, substantiated above, that the frequency 
distribution of ages at death can be derived by a linear transformation of the distribu- 
tion of vitality of the individuals at young adult age (before death sets in). Vitality was 
assumed to depend on a number of vitality variables important for survival. A single 
vitality variable may not by itself be a good “predictor” of mortality; however, if such 
a variable is sufficiently “integrated,” i.e., depends on many other physiological and 
biochemical functions, then such a variable may have a predictive value. Several con- 
ditions have to be met, however. The tested aging rate modulator should have no specific 
effect on the vitality variable, but a general effect on the body functions on which the 
variable dei>ends. It should be a relatively simple matter to measure the variable in a 
large population, to avoid exchanging “labor for time.” The method of measurement 
should be noninvasive. And, of course, it should be demonstrated hrst that in control 
populations, aging with simple passage of time affects the frequency distribution of 

vitality in a uniform and predictable way. 

It can be shown that in populations of adult mice, body weight is a vitality variable 
that meets these criteria (Economos and Miquel, 1980). Figure 23 shows the distribu- 
tion of body weight in a population of 125 male mice at various ages. A log-normal 
pattern is apparent. These mice were fed “soft” lab chow ad libitum. A population of 
45 male mice fed harder lab chow exhibited a less asymmetric distribution, with no 
overweight mice (Fig. 24); possibly, the hardness of food acted as a restraining factor 
that effected a tighter control of body weight. This group of mice was followed through 
very old age: apparently, after the age of about 24 months, the death rate begins to 
accelerate; the frequency distribution changes continuously in a poorly defined way. 
Therefore, body weight can be used as a vitality variable only in the period 8-24 
months. Studies illustrating the predictive value of body weight are documented else- 
where (Economos and Miquel, 1980). Figure 25 shot^s one example. Here, the change 
of the frequency distribution of body weight of a group of mice over a 6-month period 
during which they were fed an antioxidant-enriched diet, is compared with a change 
observed in a group of control mice. The change in the experimental group results from 
the “slimming-out” of practically all overweight mice, which is considered to be a favor- 
able effect. This correlated well with a small favorable effect of the diet on average life 
span (Economos and Miquel, 1980). 

2. Develop a method to estimate the area under the frequency distribution curve 
of a suitable vitality variable, defined in such a way that only individuals with a vitality 
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Figure 23. Age profile of body weight frequency distribution 
in a population of male mice fed “soft” lab chow ad libitum. 
From Economos and Miquel (1980). 



value above a certain threshold when tested in early adulthood should be included (Fig. 
26, top). A suitable vitality variable is, for instance, neuromuscular coordination, which 
can be estimated in simple ways both in fruit flies and in mice (Miquel et al., 1972a; 
Miquel and Blasco, 1979; Economos, 1979b). As the vitality frequency distribution 
changes with age (mainly shifting toward lower values and having larger variance), the 
percentage of individuals with vitality above the selected threshold decreases. It can be 
shown (Economos, 1979b) that if this percentage is plotted vs. age, a curve roughly 
parallel to survivorship is obtained (Fig, 26, bottom), the parallelism being closer the 
more integrated the vitality variable, i.e., the larger the number of variables on which 
it depends. Again, aging rate modulators can be expected to affect both the vitality 
curve and the survivorship curve of an experimental population in a similar way. How- 
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ever, because the vitality curve precedes the survivorship curve by up to 40% of the 
maximal life span, a time-condensing effect is achieved. Unlike in the first approach 
above, this method is suitable also during the latter part of the life span when there is 
a large death rate. 

Figure 27 illustrates the application of this approach in the study of the effect of 
two environmental factors on aging and life span of populations of fruit flies. As a vital- 
ity variable we used neuromuscular coordination. At each age, the population vitality 
(in a longitudinally studied population) was assessed based on the negative geotaxis 
drive of fruit flies (Miquel et aL, 1972a). Clearly, under both levels of each of the two 
environmental factors tested, the selected vitality variable was a good predictor of sur- 
vivorship. Finally, Fig. 28 illustrates the application of this approach to mice. On the 





Figure 24. Age profile of body weight frequency distribution in a population of male mice fed “hard” lab 
chow. From Economos and Miquel (1980). 


210 


ANGELOS C. ECONOMOS 9t «/. 


CONTROLS 
0 18 months (N = 46) 
A 24 months (N = 43) 



Figure 25. Effect of antioxidant diet on frequency distribution of body weight over a 6-month period. From 
Economos and Miquel (1980). 
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PERCENT OF POPULATION 
WITH A CERTAIN VITALITY 


Figure 26. Top: Schematic representation of the age profile of the 
frequency distribution of a vitality variable in a given population. 
[Drawn based on data on latency of first mating in fruit flics from 
Economos et al. (1979a); see also Fig. 3a.J Bottom: Percentage of 
population with vitality above a threshold value (nonhatched area 
under vitality distribution curve and above threshold line) plotted 
against age; comparison with the population's s.urvivorship curve. 
From Economos (1976b). 
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Figure 27. Population vitality curve in fruit flies as predictor of the population’s survivorship curve at two 
levels of two environmental factors. Vitality is here neuromuscular coordination assessed by the negative 
geotaxis test; the percentage of the population that pass the test is plotted vs. age. Data from Miquel et ai 
(1975, 1976). 



Figure 28. Population vitality curve in mice as predictor of the population’s survivorship curve. Body weight 
(left) and neuromuscular coordination (right) arc vitality variables considered here. (See text for more 
details.) Data from Economos and Miquel (1980) and Miquel and Blasco (1979). . 
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left, body weight was chosen as the vitality variable, and the data were from the study 
on which Fig. 25 was based. A rough parallelism was again obtained between vitality 
and survivorship curves, and the effect of the antioxidant diet on survivorship can be 
predicted from its effect on the population vitality curve. On the right, the vitality var- 
iable was neuromuscular coordination assessed by a special technique (Miquel and 
Blasco, 1979). When the data were plotted as proposed here, a vitality curve was 
obtained that again roughly paralleled the survivorship curve. 


5 . CONCLUDING REMARKS 

In this chapter we have attempted to show that “some reality lies concealed in 
biological variation” (Oldham, 1968). This “reality” has its principles, laws, mecha- 
nisms, and rules, only a few of which we have sketched. A related idea we pursued was 
that important information may be lost in the process of ignoring frequency distribu- 
tions of physiological variables (as is customary in experimental physiology and ger- 
ontology). We suggested that it may be advantageous to expand one’s “statistical field 
of vision” beyond simple averages ± standard deviations. 

Indeed, frequency distribution analysis may make visible some hidden information 
not evident from a simple qualitative analysis, particularly when the effect of some 
external factor or condition (e.g., aging, dietary chemicals) is being investigated. This 
was clearly illustrated by the application of distribution analysis in the study of varia- 
tion in mouse liver cellular and fine structure (Section 4.1), and may be true of fine 
structural studies in general 

In living systems, structure and function interact in a dynamic way; they are 
“inseparable,” unlike in technological systems or machines. Changes in fine structure 
therefore reflect changes in function. If such changes do not exceed a certain physiologic 
range, a quantitative analysis of structure will provide valuable information on quanti- 
tative changes in function that may not be possible or easy to measure directly. Because 
there is a large inherent variation in fine structure of cells in a given organ of an indi- 
vidual and among individuals, changes in fine structure can be analyzed only by study- 
ing frequency distribution curves of various structural characteristics (dimensions). 
Simple averages ± S.D. do not in general reveal all information on the effect of a 
certain factor, because often this effect is not uniform; on the contrary, this will be 
apparent from distribution analysis because the form of the curves will be affected. 

We have also attempted to show in this chapter that similar general statistical 
principles and mechanisms may be operative in biological and technological systems. 
Despite the common belief that most biological and technological characteristics of 
interest have a symmetric bell-shaped (normal or Gaussian) distribution, we have 
shown that more often than not, distributions tend to be asymmetric and often resemble 
a so-called log-normal distribution. We saw that at least three general mechanisms may 
be operative, i.e., nonadditivity of influencing factors, competition among individuals 
for a common resource, and existence of an “optimum” value for a studied character- 
istic; more such mechanisms could exist. As we stressed, all these mechanisms are 
purely statistical and thus very general, which is illustrated by the similarity of fre- 
quency distribution curves of structural and fine structural characteristics of cells with, 
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e.g., the curves for physiological vitality, mortality, body weight, latency of mating of 
flies, glomerular basement membrane width of hamster kidneys, age of fellows of the 
Royal Society of London, professorial salaries, and capacitor breakdown voltage, to 
mention only some of the examples we have discussed. 

In conclusion, biological variation should be fully recognized and accepted and the 
effect of aging thereupon should be studied. Thus, variation may become an investiga- 
tive tool for the experimental gerontologist: when appropriately studied, it can speed up 
and thus increase the scope of much of experimental gerontology. 
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